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Abstract 

Networks representing many complex systems in nature and society share some common structural properties 
like heterogeneous degree distributions and strong clustering. Recent research on network geometry has 
shown that those real networks can be adequately modeled as random geometric graphs in hyperbolic spaces. 
In this paper, we present a computer program to generate such graphs. Besides real-world-like networks, the 
program can generate random graphs from other well-known graph ensembles, such as the soft configuration 
model, random geometric graphs on a circle, or Erdos-Renyi random graphs. The simulations show a good 
match between the expected values of different network structural properties and the corresponding empirical 
values measured in generated graphs, confirming the accurate behavior of the program. 
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Program summary 

Program title: Hyperbolic Graph Generator 

Program summary URL: http: // named-data. github. 

io/ Hyperbolic-Graph-Generator/ 

Licensing provisions: General Public License, version 3 
Programming language: C+ + 

Computer/Operating system: Any 

Nature of problem: Generation of graphs in hyperbolic 
spaces. 

Solution method: Implementation based on analytical 
equations. 

Running time: Depends on the number of nodes. A few 
seconds for graph in the example provided. 

Other features: Can be used as a command-line tool or 
installed as a library to support more complex software. 
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chiara@caida.org (Chiara Orsini), dima@neu.edu (Dmitri 
Krioukov) 


1. Introduction 

The interactions between components of a complex 
system are often represented as a network. This mod¬ 
eling allows for rigorous mathematical treatment, 
and broadens our understanding of the system [1]. 
Many real networks possess common structural pat¬ 
terns, including, in the first place, heterogeneous 
(often power-law) distributions of node degrees [2], 
and strong clustering, i.e., higher numbers of tri¬ 
angular subgraphs than predicted by classical ran¬ 
dom graph models [3]. Recently introduced geo¬ 
metric graph models, based on the assumption that 
nodes in real networks are embedded in latent hyper¬ 
bolic spaces [4, 5], reproduce these common struc¬ 
tural properties of real networks. Furthermore, these 
hyperbolic graphs replicate dynamical processes on 
top of real networks [6] and accurately predict miss¬ 
ing links in them [7]. 

In this work we present a program to generate ran¬ 
dom hyperbolic graphs. This software implements 
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Table 1: Regimes in the model. 


and extends the network model introduced in [4]. 
Nodes are randomly sprinkled on a hyperbolic disk, 
and the probability of the existence of an edge (the 
connection probability) between two nodes is a func¬ 
tion of the distance between the nodes in the hy¬ 
perbolic space. Thus generated graphs have strong 
clustering, and node degree distributions in them 
are power laws. Moreover, other popular and well- 
studied random graph ensembles, namely the soft 
configuration model (SCM) [ 8 ], (soft) random geo¬ 
metric graphs (RGGs) on a circle [9, 10], and Erdos- 
Renyi (ER) random graphs [11], appear as degenerate 
regimes in the model. Table 1 shows all the model 
regimes, the total of six. Each regime is defined by 
the values of only two parameters: 7 , which is the 
expected exponent of the power-law degree distribu¬ 
tion, and temperature T, the parameter controlling 
the strength of clustering in the network. 

Researchers in different disciplines may benefit 
from the use of random hyperbolic graphs in their 
work. Yet the full implementation of the model and 
all its regimes is a tricky business, which involves 
dealing with some delicate details, due to a variety 
of internal parameters and their interactions over the 
six regimes. In Sections 2-4 we describe the imple¬ 
mentation details of the model, including how all the 
parameters are calculated in each regime. A good 
match between the values of expected graph proper¬ 
ties and their observed values in generated graphs is 
confirmed in Section 5. 


2. Input parameters and coordinates 


combination of 7 and T values will define the graph 
ensemble from which generated networks are sampled 
(Table 1). 

Given the input parameters, the graph generation 
process consists of three steps: 

1. Compute the internal parameters, such as the ra¬ 
dius R of the hyperbolic disk occupied by nodes, 
as functions of the input parameter values, Sec¬ 
tions 3, 4. 

2. Assign to all nodes their angular and radial co¬ 
ordinates on the hyperbolic plane, Section 2.2. 

3. Connect each node pair by an edge with prob¬ 
ability (the connection probability), which is a 
function of the coordinates of the two nodes, Sec¬ 
tions 3, 4. 


2.2. Coordinate sampling 

The assignment of node coordinates is done as fol¬ 
lows in all the six regimes. 

Angular coordinates 9 of nodes are assigned by 
sampling them uniformly at random from interval 
[0, 27 t), i.e., the angular node density is uniform 
p{6) = 1/(2?r). 

Radial coordinates r G [0, R], where R is the radius 
of the hyperbolic disk, are sampled from the following 
distribution, which is nearly exponential with expo¬ 
nent a > 0 , 


p{r) = 


sinh( 


cosh aR — 1 


x(r—R) 


(i) 


2.1. Input parameters 

The program input parameters are the number of 
nodes N, the target expected average degree k of the 
network, the target expected power-law exponent 7 
of the degree distribution, and temperature T. The 


The calculation of internal parameter R is described 
in detail below; it is different in different regimes. 
Internal parameter a depends on the expected ex¬ 
ponent 7 of the power distribution P(k ) ~ fc -7 of 
nodes degrees k in generated graphs, and on the cur¬ 
vature of the hyperbolic space £ = \J—K. which is 
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set to ( = 1 by default. For temperatures T < 1, this 
relationship is given by 

7 = 2 “+l, ( 2 ) 

while for T > 1 it becomes 

7 = 2“T + 1. (3) 

To sample radial coordinates r according to the dis¬ 
tribution in Eq. (1), the inverse transform sampling 
is used: first a random value U t is sampled from the 
uniform distribution on [ 0 , 1 ], and then the radial co¬ 
ordinate of node i is set to 

Vi = — acosh (1 + (coshai? — 1 ) UA , for i = 1, N. 
a 

( 4 ) 

3. Regimes with finite 7 > 2 

3.1. T € (0,oo): Soft hyperbolic random geometric 
graphs 

This is the most general regime in the model, from 
which all other regimes can be obtained as limit cases. 
The connection probability in this case is 

P^= i + e W72)(*-fl) ’ {5) 

where /3 = 1/T, and R is the radius of the hyper¬ 
bolic disk occupied by nodes. The hyperbolic dis¬ 
tance x between two nodes at polar coordinates (r, 9) 
and (r', 0 ') is given by 

x = ^arccosh (cosh (r cosh fr' 

— sinh sinh Cr 7 cos AO) , (6) 

where A 9 = 7 r — | 7 r — \9 — 0'\\ is the angular dis¬ 
tance between the nodes. To calculate the expected 
degree of a node at radial coordinate r, without loss 
of generality its angular coordinate can be set to zero, 

9 = 0, so that its expected degree can be written as 

N r R r 

k(r) = — / p(r') / p{x) dO'dr'. (7) 

n Jo Jo 


The expected average degree in the network is then 

- f R 

k= p(r) k{r)dr 

R R ( 8 ) 
N r R r R r 

= — / p(r) / p(r') / p{x) d0'dr dr. 

n Jo Jo Jo 

Given user-specified values of input parameters N , 
fi = 1/T, ( and k , the last equation is solved for R 
using the bisection method in combination with nu¬ 
meric evaluation of the integrals in the equation. The 
MISER Monte Carlo algorithm from the GSL library 
is used to compute the multidimensional integral in 
Eq. ( 8 ). The iterative bisection procedure to find R 
stops when the difference between the value of the 
computed integral in Eq. ( 8 ) and the target value of 
k is smaller than a threshold that is set to 1 (R 2 by 
default. 

3.2. Limit T —» 0: Hyperbolic random geometric 
graphs 

In the T —> 0 (/3 —> 00 ) limit, the connection prob¬ 
ability in Eq. (5) becomes 

p{x) = 0(R- x), (9) 

where @(x) is the Heaviside step function, meaning 
that two nodes are connected if the hyperbolic dis¬ 
tance x between them is less than i?, or they are not 
connected otherwise. The expected average degree in 
the network is given by the same Eq. ( 8 ), but with 
p(x) in the last equation. The value of R is deter¬ 
mined using the same procedure as in Section 3.1. 
The only difference is that function p{x) is given by 
Eq. (9). 

3.3. Limit T —> 00 : Soft configuration model 

According to Eq. (3), in the T —> 00 limit with 
finite a, to have finite 7 , curvature should also go to 
infinity, £ —> 00 , such that rj = £/T is finite, and 
instead of Eq. (3) one gets 

7 = 2— + 1. (10) 
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More importantly, one can show that as a result 
of £ —> oo, the expression for hyperbolic distance x 
between two nodes in Eq. ( 6 ) degenerates to 

x = r + r', ( 11 ) 

meaning that in the T —> oo regime the angular coor¬ 
dinates are completely ignored. The connection prob¬ 
ability becomes 

P( r , r ') = i + e (7)/2)(r+r'—R) ’ (12) 

and the expected average degree in the network is 

/*R rR 

k = N p(r) / p(r') p(r,r') dr'dr. (13) 

Jo Jo 

The value of R is determined using the same combi¬ 
nation of the bisection method and numeric integra¬ 
tion as in the previous section, except it is applied to 
Eq. (13). 


4. Regimes with infinite 7 —7 00 

While in the T —>• 00 limit the angular coordinates 
are ignored, in the 7 —> 00 limit the radial coordi¬ 
nates are ignored. One can show this formally by ob¬ 
serving that in this limit the radial node density ap¬ 
proaches a delta function -all nodes are placed at the 
boundary at infinity of the hyperbolic plane, mean¬ 
ing that only angular coordinates determine distances 
between nodes. 

4-1- T £ (0, 00 ): Soft spherical random geometric 
graphs 

In this most general case with infinite 7 , one can 
show that the connection probability in Eq. (5) de¬ 
generates to 


p(M') 


1 


(14) 


where Ad = 7 r — \to — \9 — 0'\ \ is the angular distance 
between the two nodes as before, while A is a param¬ 
eter controlling the average degree k in the network, 


analogous to R in the regimes with finite 7 . With¬ 
out loss of generality we can set 9 = 0 , so that the 
expression for k is 

-NT 1 

k = - - dd' =N 2 F 1 (1 1 T-T+ 1;-A), 

* •'0 i + A ( e ~y 

(15) 

where 2 Ei is the Gauss hypergeometric function, and 
T = 1//J. In the special case with T = 1, the last 
expression simplifies to 


k _ log(l -I- A) 
N ~ A 


(16) 


If T 7 ^ 1, the hypergeometric function in Eq. (15) 
cannot be evaluated using the GSL library, because 
the 2 Ei evaluation in the library is implemented only 
for the case where the fourth argument of the func¬ 
tion (—A in Eq. (15)) is between —1 and 1, while for 
sufficiently large N/k, A is always larger than 1 in 
Eq. (15). To avoid this difficulty, the following trans¬ 
formation is used [ 12 ]: 


- = 2 *i(l,T;T+l;-A) 


T 


A+1T-1 


2 F\ (1,1; 2 — X 1 ; 


1 )+ 1 


7 tT 


A + 1 


A T sin 7 tT ’ 
(17) 


If T > 1 is an integer, the second term in (17) diverges 
due to the sin function in the denominator, while the 
first term diverges because the third parameter of 
the 2 Ei function is a non-positive integer. Hence, for 
integer values of temperature T > 1, their value is 
approximated by T + e, where e is set to 10 -6 by 
default. The error caused by this approximation is 
negligible. Equation (17) (or (16) if T = 1) is then 
numerically solved for A using the bisection method, 
yielding the target value of k in Eq. (15). 


4-2. Limit T —» 0: Spherical random geometric 
graphs 

One can see from Eq. (17) that the solution for A 
at small T < 1 scales with N/k as A = (N/k)@, /3 = 
1/T. Therefore for /3 1 the connection probability 

in Eq. (14) can be written as 


p(M') 


1 

1 + (i¥f' 


4 


(18) 
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(a) Average degree. 
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Table 2: Observed properties in generated graphs (mean ± std.dev for 10 3 graph samples) with target 
average degree k = 10 . 


which in the /3 —> oo limit becomes 


5. Simulations 


p(M') = e 



NA6\ 

kir ) ’ 


(19) 


meaning that two nodes are connected if the angular 
distance A 9 between them is smaller than nk/N, 

A9<tt^, ( 20 ) 

or they are not connected otherwise. This connectiv¬ 
ity threshold ensures that the expected average de¬ 
gree in the network is k. 


4-3. Limit T —> oo: Erdos-Renyi graphs 

In this most degenerate regime, both angular and 
radial coordinates are completely ignored. This 
regime is formally achieved by keeping both a and £ 
finite while letting T —> oo. One can then show that 
the connection probability in Eq. (5) degenerates to 

P( x ) = TT^y ’ ( 21 ) 

which for sparse graphs with k <C IV tends to 
p(x) = k/N, i.e., the connection probability in clas¬ 
sical (Erdos-Renyi) random graphs. 


Tables 2a and 2b show, respectively, the average 
degree and clustering values in generated graphs for 
different regimes, 10 3 samples in each regime. All 
the regimes match the target average degree k = 10 , 
although low values of 7 lead to much higher fluctu¬ 
ations. For any 7 > 2, average clustering decreases 
with temperature from a maximum at T = 0 to zero 
at T —» 00 . 

Figure 1 shows the observed distribution of node 
degrees for three different values of 7 . As expected, 
for finite 7 , distributions follow a power-law P{k) ~ 
A; -7 . In the limit at 7 —► 00 the observed degree 
distributions follow exactly the Poisson distribution 
with mean p = k. 

Figure 2 shows clustering in generated graphs. For 
finite 7 , low-degree nodes have stronger clustering 
than high-degree nodes. In the case of 7 —> 00 , all 
nodes have similar clustering values. 
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Figure 1: Observed degree distributions in generated 
networks of size N = 10 4 and target average degree 
k = 10. The results are averaged across 10 3 generated 
graphs for each combinations of the parameters. 


Figure 2: Observed clustering in generated networks 
of size N = 10 4 and target average degree k = 10. 
The results are averaged across 10 3 generated graphs 
for each combinations of the parameters. 
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